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Abstract 

We propose a traffic flow model in which the vehicles are filed from their maximal ve- 
locities, the fast cars run with Vmaxl, whereas the slow ones run with Vmax2. Using new 
overtaking rules which deals with deterministic NaSch model, it is found that the funda- 
mental diagram exhibits three new topologies, depending on the fractions ff as t an d fslow 
of fast and slow vehicles respectively, in which the current profile displays two branches 
with negative slopes and two branches with positive ones. Moreover, in the second branch 
of the fundamental diagram, the model exhibits an absorbing phase transitions in which 
the behaviour of the order parameter and the current J is described by the power laws. 
In this case, it 'is found that the system present a universel scaling law. On the other 
hand, a simple change in the rule of overtaking induce the metastability which depends 
on the state of the chain instead of external parameters [4, 5, 6]. Furthermore, in the case 
of random fractions of vehicles, the fundamental diagrams are similar to the experiments 
results [7, 8] 
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1 Introduction 



The main raison for studying vehicular traffic is its no stop increasing problems. Indeed, a 
new study suggests that traffic fumes is killing tens of thousands of people across Europe, 
and the costs of treating the illness associated with traffic pollution exceeds the costs 
arising from traffic accidents]!]. On the other hand, the congestion traffic has a harmful 
impact on the economy since drivers and goods waste more time and fuel standing in 
a traffic jams( The Texas transportation institute found that the national cost of traffic 
congestion is 62.3 billion $ annually) [2] . Moreover, the construction of new highway is 
not affordable and publicly acceptable due to the expensive costs and for want of lands 
especially in urban areas. Many solutions and strategies were proposed in order to use 
the available infrastructures more efficiently, but with little success. For such raisons 
scientists including physicists, mathematicians, chemists and engineers are still working 
in this field in order to give a general solution, which can reduce the damage caused by 
roadway traffic. 

Obviously , traffic flow depends strongly on the density C of cars on the road, but how 
vary the flux J and the average velocity V with vehicles density C? Functional relations 
between J, V, and C have been measured since 1935 by Greenshields (see[9]). However, it 
is difficult to obtain very reliable and reproducible detailed empirical data on real traffic 
for several reasons namely: It is not possible to perform such laboratory experiments on 
vehicular traffic, thus empirical data are collected through passive observations rather than 
active experiments. Secondly the interpretation and evaluation of the collected data is not 
a simple task because traffic states depend on several influences like: weather conditions, 
presence of incidents, the measurement conditions ( location where the traffic variable 
were measured, speed limits, ...), and other irregularities [9, 10] . As a consequence, the 
empirical determination of the dynamical properties of traffic flow is not evident. In order 
to use the empirical results for a theoretical analysis, many authors have used the mean 
values of the flow at given density which lead to a collection of possible forms of averaged 
fundamental diagrams(J,C) [3, 7] consistent with empirical data. The availability of these 
data has stimulated a great number of experimental studies and their comparison with 
theoretical models. The physical theory of traffic is based on the concept and techniques 
of statistical mechanics [11, 12, 13], the complex behaviour of traffic flow on a freeway has 
been the major topic of current interests [14, 15]. Indeed, many theoretical research using 
both analytical and numerical simulations in one and two dimension, has been done in 
order to make clear the characteristics of this complex behaviour [16, 17]. In this context 
there are two different ways for modelling traffic: The macroscopic models which are 
based on fluid-dynamical description and the microscopic ones where attention is explicitly 
focused on individual vehicles which are represented by particles. The interaction is 
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determined by the way the vehicles influence each others movement. In other words 
in the microscopic theories traffic flow is considered as a system of interacting particles 
driven far from equilibrium [9]. Thus, it offers the possibility to study various fundamental 
aspects of the dynamics of truly non equilibrium systems which are of current interest 
in statistical physics[9, 18, 19]. Within the conceptual framework of the microscopic 
approach, the particle hopping models describe traffic in terms of stochastic dynamics of 
individual vehicles which are usually formulated using the language of cellular automata 
(CA)[20]. In general, CA are an idealization of physical systems in which both space 
and time are assumed to be discrete and each of the interacting cell can have only a 
finite number of discrete states, thus in CA models of vehicular traffic first proposed by 
Nagel and Schreckenberg[21]and subsequently studied by other authors using a variety of 
techniques [22, 23], the position, speed, acceleration as well as time are treated as discrete 
variables and the lane is represented by a one-dimensional lattice, each site represents a 
cell which can be either empty or occupied by at most one vehicle at a given instant of 
time. Hence, these properties lead to a complex dynamic behaviour and a fast simulation 
of a great number of interacting vehicles. 

However, the vehicles on the road run with different speed. Thus, for modelling traffic with 
different kind of vehicles, for example cars and trucks, we consider ff ast the fraction of fast 
vehicles(typel), i.e, vehicles with a higher maximal velocity V max i, while the remaining 
fraction I — ffast is assigned to slow vehicles(type2), i.e, vehicles with low maximal velocity 
V ma x2- Hence, for a fixed density C, we denote by C fast = C f fast and C s i ow = C f siow 
the density of typel and type2 respectively. Such situation leads to the formation of 
coherent moving blocks of vehicles each of which is led by a vehicle with low V max , i.e., 
the phenomena of platoon formation [24] . The strong effect of slow vehicles is also shown 
in the multi-lane systems [25, 26], indeed, the simulation results revealed that already 
for small densities, the fast vehicles move with the average free flow velocity of the slow 
cars even if only a small fraction of slow vehicles have been considered. This effect is 
quite robust for different choices of CA model as well as for different lane changing rules. 
In this paper, we propose a CA traffic flow model which we hope can contribute to a 
better understanding of the traffic flow. Our analysis is inspired by the situation in real 
traffic hence, the proposed model is based on the following facts: a mixture of fast and 
slow vehicles taking into account the disorder in their maximal velocities, more real traffic 
overtaking rules which deals with deterministic NaSch model in a single lane[21, 23], 
met ast ability, and random fractions of different kind of vehicles. For didactical reasons 
the paper is organized as follow; section 2 is devoted to explain the model. In section 3 we 
will present the main results obtained by simulation with a critical discussions including 
an analytic treatment, scaling analysis and the metastability. In the last section we will 
present a general conclusion. 
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2 Model 



It is well known that Nagel and Schreckenberg are among the early pioneers who have 
addressed the problems of traffic flow; their CA model[21] has stimulated a great numbers 
of publications (for review see[9, 10]). Indeed, the NaSch model is based on four steps 
which are necessary to reproduce the basic features of real traffic; however more addi- 
tional rules should be added in order to capture more complex situations. In this context 
we have formulated our model. Indeed, the case of disorder on the maximal velocities 
leads to platoon phenomena as was discussed in the introduction, this means that the fast 
vehicles(typel running with V max i) are trapped in a coherent moving blocks which are led 
by a slow vehicle (type2 running with V rnaX 2). Thus, in order to overcome this situation 
fast vehicles must overtake the slow ones, which will lead to a competition between the 
formation of platoons and their dissociation by overtaking. However, the choice of the 
overtaking rules is not a simple task due to the strong effect of slow vehicles. The situation 
of overtaking is considered only when the following rules (which deal with deterministic 
NaSch model) are satisfied: 
CI: typel follows type2 

C2: Gapl < Vmax2 

Gapl is the distance between the two vehicles typel and type2. This choice ensures that 
fast vehicles driving in a slow platoon try to overtake if it is possible. 

C3: Gap2 > gap c 

Gap2 is the number of empty sites in front of type2, i.e., the gap between type2 and 
the next vehicle in the chain, while gape = min(V \type2) + 1, Vmax2) + 1 is the min- 
imal safety distance required for overtaking. Indeed, the vehicles in the chain move at 
the same time step(parallel dynamic), and according to the deterministic NaSch rules 
V(type2) = min(V(type2) + l,Vmax2,Gap2), hence, for large Gape2 the velocities of 
type2 in the next time step will be V(type2) = min(V(type2) + l,Vmax2), as a conse- 
quence C3 ensures that at the situation of overtaking, type2 moves without any hindrance 
which means that typel estimates the velocity of type2 before overtaking. 

C4: if p s > rand, then typel overtakes with V {typel) = min(NGap, Vmaxl). 
p s is the overtaking probability and Rand is a random number. On the other hand, to 
avoid collision with vehicles in front of type2: NGap=Gapl+l+Gap2 is the new gap re- 
quired for overtaking which is obtained according to C2 and C3. Thus when the above 
criteria are satisfied typel overtakes with its maximal velocity if enough empty space is 
allowed. 
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C5: Vmaxl > 2{Vmax2 + 1) 

This inequality guides the choice of Vmaxl and Vmax2; it is obtained according to the 
following criteria C2, C3 and C4. 

In general, the update in this model is divided into two sub steps: 

in the first sub step (I), typel which satisfies the four criteria can overtake in parallel; 
However, in the second substep (II), the remaining vehicles (i.e, typesl which didn't sat- 
isfy the four criteria and types2) move forward in parallel according to the NaSch rules, 
without taking into account the new positions of the typel moved in the first sub step. 

3 results and discussion 
3.1 Simulation results 

In our computational studies we have considered a chain with N = 1000 sites and periodic 
boundaries. In this case simulation begins with particles randomly distributed around the 
chain according to their densities Cf ast = Cff ast and C s i ow = Cf s \ om where C is the global 
density. The systems run for more than 40 000 MCS to ensure that steady state is reached, 
at this moment data including the current and velocity are collected. In order to eliminate 
the fluctuations about 100 initial configurations were randomly chosen. 

For p s = 1.0 i.e., deterministic case and depending on the values of ff as ti the fun- 
damental diagram (J,C) exhibits three new topologies in which four different branches 
occur. In fact, at free flow we have two branches: one with positive slope and the other 
with negative one, which meet at local maximal current. The same behaviour occurs at 
congested traffic flow. The maximal current at free flow can be smaller, equal or greater 
than the one located at congested traffic flow. Indeed, for Vmaxl = 5 and Vmax2 = 1 
in Fig.l-a the fraction of the fast vehicles is smaller than the fraction of the slower ones 
ffast « i slow 1 then the maximum of the free flow branch is smaller than the one of the 
congested branch. For ff as t >> fsiowi the maximum of the free flow branch becomes 
greater than the congested branch (Fig.l-c), while for ffast ~ fslow the two maximums 
are equal (Fig.l-b). Such topologies in the fundamental diagram was obtained in em- 
pirical work[3]. Moreover, in Figs.l, we can distinguish several different states of traffic 
depending on the value of the density C of vehicles in the road: 

For C < C max \ (the first branch), the current increases with C to reach at C max \ its 
first maximal value J m axi which depend strongly on the values of ffast , indeed, {C max \ , J ma xi ) = 
{(0.14,0.288); (0.16,0.48); (0.19,0.68)} for f fast = {0.3,0.6,0.75} respectively. In this re- 
gion, the four criteria of overtaking are satisfied due to the lower values of C, i.e, an 
enough empty space between vehicles, thus all vehicles typel overtake type2 which means 
that both types of vehicles are in their free flow regime as shown in Fig. 2. 
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For C max \ <C< C m i n (the second branch); C m in is the density at which the current 
show a nonzero minimum J m in- in this region the current value fall with C. Indeed, the 
situation of overtaking becomes more difficult as long as C increases, thus not all vehicles 
typel can overtake, which leads to a competition between the formation of platoons and 
their dissociation by overtaking; as a consequence typel are in their congested phase 
whereas type2 are still in their free flow regime as displayed in Fig. 2. 

For C m i n <C< C maX 2 (the third branch) typel can not overtakes type2, such situa- 
tion leads to the formation of coherent moving blocks of vehicles each of which is led by 
type2, i.e, the phenomena of platoon formation. 

We recall that C maX 2 is the density at which the current show a second maximal value 
Jmax2- Thus for C > C maX 2{th.e fourth branch) both vehicles are in the jamming phase. 

In Fig.2 we have plotted the mean speed of each types of vehicles versus the density 
of cars C with the same fractions as in Fig.l-c. Hence Fig.2 proves that for (C < C max i) 
(the first branch) each types of vehicles run with their average free flow velocity. For 
Cmaxi <C< Cmin (the second branch) the average speed of typel fall with C while the 
average speed of type2 remain constant, which means that types2 are still in their free 
flow regime whereas typesl are in their congested phase. For C m in < C < C ma x2 (the 
third branch) typel and type2 run with the same speed .i.e., typel drive in a slow platoon 
which is led by type2; While for C > C ma x2 the jamming phase take place. 

We conclude that for typel the first maximum (C max i, J ma xi) corresponds to the 
transition from free flow to congested phase. (C m j n , J m in) indicates the transition from 
congested phase to platoon phase while the second maximum (C ma x2, Jmax2) corresponds 
to the transition from platoon phase to jamming phase of both types of vehicles. 
On the other hand, the same analysis hold for the case when p s < 1.0 except that J m axi 
is lowered whereas C ma xi and C m in are shifted to the left as displayed in Fig.3. Indeed, 
even if enough empty space exist not all typel can over take due to C4. 
Moreover, J m in does not depend on ff as t and p s , in contrast to J ma xi- This is due to the 
structure of the chain in this region, i.e., there is no sufficiently empty space for overtaking 
in comparison to the first region. Also C > C m in imply the end of the competition, hence, 
the movement of typel is dictated by type2 which leads to the well known fundamental 
diagram of one species of vehicles with V max = 1, for such raison J maX 2 is also independent 
on ff^ and p s . 

Furthermore we believe that our model can be extended to the case of n types of vehicles 
with n different maximal velocities. In this context we have computed the case of three 
types of vehicles, where f fasti, ffast2 and f siow are the fractions of typel (V max i = 10), 
type2 (V maX 2 = 4) and type3 (V maX 3 = 1) respectively. The choice of V max follows the 
condition C5 mentioned above. In this case typel can overtakes type2 and type3 whereas 
type2 can overtakes only type3. Hence, the current displays three local maximums. At 



6 



free flow the maximums are comparable because both typel and type2 are fast. Indeed, 
Fig. 4 shows that for low density the three vehicles are in their free flow regime, however, 
the mean speed of typel falls quickly with C in comparison with the average speed of 
type2, which explains the appearance of the first and second maximum. As long as C 
increases typel and type2 are not able to overtake, hence, they drive in a slow platoon 
which is led by type3. 



3.2 Analytic treatment 

As was discussed in the previous section the disorder in the maximal velocities and the 
overtaking rules lead to a competition between the platoons formation and their dissoci- 
ation. Thus we need a physical quantity to quantify this competition. 
Lets consider Nf the number of vehicles typel in the steady state, for which the rules 
C1,C2 and C3 are satisfied .i.e., fast vehicles which are favourable to overtake. Unfortu- 
nately, not all the Nf vehicles can do because C4 will not be fulfilled for all these cars. 
As a consequence, among Nf only I vehicles succeed to overtake depending on the value 
of p s . The probability of such event is given by: 



P=C I NfP I s (l-p s ) N f- 1 (1) 



With 



C 



j N f ! 
N f (N f -I)!I! 



Hence, the mean value of I is: 

<i>=j2ic i NfP i(i- Ps ) N f- i (2) 



v 



However, 
Thus, 



i(? N =N f CZU (3) 



N f 



< I >= N fPs ]T C^plHl - Ps )W-i)-(*-i) ( 4 ) 
1=1 

for i — > i — 1, the sum becomes a binomial law: 

N f -1 

E Cjr / -irt(l-P-) (JV/ - 1)_i = l (5) 



As a consequence, 

<I>=N fPs (6) 

On the other hand, we denote by fd the fraction of vehicles typel which succeed to 
overtake: 

N f 

fd = T r Ps (7) 
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Where, N p is the number of vehicles typel which are trapped behind a slow vehicles 
.i.e., the number of fast vehicles for which at least CI is satisfied. Thereafter, f d is 
the physical quantity which quantify the competition between overtaking and platoons 
formation, indeed: 

For C < C m axi (the first branch)where there is enough empty space for overtaking, fd 
vehicles typel run with V ma xi whereas the remaining fast vehicles (1-fa) moves with the 
same velocity V maX 2 as type2; Thus the current J can be expressed as: 

J(C) = Cf ast [fd V m axl + (1 — fd) Vmax2\ + C s l ow V max 2 (8) 

We recall that Cf as t = C ff as t and C s i ow = C f s i ow are the density of typel and type2 
respectively, and ff as t = 1 — fsiow, hence, one can easily get: 

J(C) = C [ff as t fd(Vmaxl ~ V ma x2) + V max 2] (9) 
Furthermore, in this range of density, Nf = N p which leads to 

fd = Ps (10) 

Hence, the current in the first branch is: 

J(C) = C [ff as t Ps{Vmaxl ~ Vmax2) + Vmax2\ (H) 

For C m in < C < C ma x2 (the third branch) where overtaking is note allowed, Nf « 1 
thus 

fd = (12) 
which leads to current in the platoon phase: 

J(C) = C V maX 2 (13) 
For C > C maX 2(the fourth branch), here the current is limited by the density of holes: 

J(C) = 1-C (14) 

which is a characteristic of the jamming phase. However it is difficult to do analytical 
calculation for NfandN p in the second branch C ma xi < C < C m i n , indeed the simulation 
results Fig. 5 shows that fd is a singular function in this range of density even in the 
deterministic case p s = 1., hence, the current in this branch keep the formula (9) where 
fd can be obtained by simulation. Furthermore, we remark that the expressions of the 
current above are obtained depending on the value of fd in each branch except the last 
one (jamming phase), thereafter, the fundamental diagram(J,C) Fig.l-c is given by: 

J(C) = Mifl{C[f d ffast{V m axl - V max2 ) + V max2 ], 1 - C} (15) 
This general formula sum up the previous expressions, however, it'is not possible to obtain 
(Cmaxi, Jmaxi) and (Cmin, J min) by a simple derivative due to the fact that fd is a singular 
function in the second branch. Thus, in the following section we will present a new way 
to study this singular behaviour. 
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3.3 Power laws and scaling analysis 

It' is well known that in equilibrium systems, there is a certain lower dimension below 
which the system is always disordered. However, the situation is very different in nonequi- 
librium systems. This fact has initiated intense research regarding the phase transitions 
in one-dimensional driven diffusive systems far from equilibrium [9]. Moreover, the con- 
cepts and techniques developed for equilibrium phase transitions can also be applied to 
nonequilibrium systems such as the order parameter concept, scaling, universality as well 
as renormalization group, therefore, numerous nonequilibrium critical phenomena were 
successfully investigated by such methods [27]. Here we are interested in the so called Ab- 
sorbing Phase Transitions(APT)[27, 28, 29, 30] which occur in nonequilibrium dynamical 
systems that are characterized by at least one absorbing state. Indeed, in our model the 
platoon phase(the third branch) is an absorbing state in which the fast vehicles become 
trapped forever, whereas the active phase correspond to free flow phase of typel(the first 
branch). The competition between the proliferation and the annihilation process which 
constitute the essential physics of the AFT, is described in this model by the competi- 
tion between the formation of platoons (proliferation process) and their dissociation by 
overtaking(annihilation process). However, in the active phase(C < C max \) the annihila- 
tion (dissociation of platoons) prevail the proliferation(platoons formation) thus, fd the 
physical quantity which quantify this competition have a non zero v&lue(fd = p s )- In 
the absorbing state(C m j n < C < C max 2) in which proliferation outweight annihilation, 
fd = 0. For C max i < C < C m i n the system undergoes an absorbing phase transition from 
the active phase to the absorbing phase where fd fall with a singular manner as long as C 
approach the critical point C m i n . As a consequence, fd represent the order parameter of 
the system and the density C is the control parameter of the absorbing phase transition. 
Moreover, it' is well known that critical systems are characterized by power laws which 
describe the singularities observed sufficiently close to the critical point; fortunately the 
singularity in this model appears in a specific region of the density i.e., when C max \ < C < 
C m i n , furthermore, Fig.5 show that this region depend only on ff as t] Thus the singular 
part of the order parameter is assumed to scale asymptotically as: 

fd(C) ~ f] ast g(C/f fast ) (16) 

Where g(x) = cte is a scaling function, which means that the singular part of fd is 
asymptotically a generalized homogeneous functional]. Indeed, the simulation results 
Fig.6 confirm this assumption: 
For f fast > 0.6 

7 = and 5 = 0.773 

However, for ff as t < 0.6 the singular part of fd is independent of ffast as shown in inset 
of Fig.6. i.e., both 7 = and 5 = 0. 
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On the other hand, the same analysis hold for the second branch of the fundamental 
diagram. Indeed, in section 3.1 we have shown that (C max i, Jmaxi) depend on both 
ffast an d p s , as a consequence the current in the second branch is also asymptotically a 
generalized homogeneous function, hence for p s = 1 it' is assumed to scale asymptotically 
as: 

J(c) ~ ff ast f(c/r fast ) (17) 

with the scaling function 

f(x) = cte 



This result is confirmed by simulation as displayed by Fig.7. However, the exponent 
a and v depend strongly on the value of ff as t- 

For ff a st > 0.6,i.e., when the vehicles type2 are considered as an impurity ffast » 

f slow- 
est = 1.631 and u = 0.773 

For 0.4 < ffast < 0.6, here both type of vehicle are comparable ffast ~ fslow- 

a = 1.3 and v = 

For ffast < 0.4, in this case the vehicles typel are now considered as an impurity 

ffast fslow- 

a = 0.532 and v = 

However, for a fixed fraction of typel ff as t an d p s < 1 (the stochastic case) Fig.7, the 
current scale asymptotically as: 

(8) J(C) ~ p z s h{C/p1) (18) 

with the scaling function 

h(x) = cte 



Here we can distinguish two regions depending on the value of p s , indeed: 
For 0.8 < p s < 1 (p s — > 1) i.e., the system tend to the deterministic case: 

z = 0.9613 and n = 1.0591 

Whereas for p s < 0.8: 

z = 0.302 and n = 

According to these results we can formulate some remarks: 

Firstly, the critical exponents 7, 5, a, u, z, n are not universal because they depend on 
the parameters of the systems. However, the model shows a universel scaling law S = v , 
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which confirm that the behaviour of the system in the second branch of the fundamental 
diagram is imposed by the one of fd in this range of density. Secondly, ff as t = {0.6; 0.4} 
represent a crossover, indeed, the system change its critical behaviour depending on which 
types of vehicles is the majority, such result was mentioned in beginning of section 3.1. 
Furthermore, the strong effect of p s appears only when p s < 0.8. 

3.4 Metastability 

According to the previous sections, there is a density region C max \ < C < C m j n (the 
second branch) in which we can have both the free flow and the congested traffic Fig.2. 
Such remark indicate that the system can exhibits metastability and the related hysteresis 
effects [7]. This phenomena have been observed in coupled-map lattice models [23], and 
in some generalizations of the NaSch model which are based on slow to start rules[4, 5, 6]. 
However, in the latter cases, the slow to start rule is formulated by a modifications of 
the braking steps of the original NaSch model, whereas in our model it is obtained by a 
simple change in the second sub step (II) of updating. Indeed, when C3 is not satisfied, 
the corresponding vehicle typel stop in this time step, it will wait until the next time 
steps to move forward according to NaSch rules. Hence our slow to start rule depends on 
C3 ,i.e, on the configuration of the chain instead of external parameters as in [4, 5, 6]. 
Starting from two different initial conditions, the mega jam (a large compact cluster of 
standing vehicles) and the homogeneous state (vehicles are distributed periodically with 
the same constant gap) we obtain the hysteresis observed in Fig.8. 

4 Conclusion 

The investigations made in the above sections lead to the following conclusion: 
We have studied the case of mixture of fast and slow vehicles taking into account the 
disorder in their maximal velocities, the key ingredient of our proposed model is how the 
faster cars can overtake the slow ones. Indeed, using a new overtaking rules which deal 
with deterministic NaSch model in a single lane, we have found that the proposed model 
belong to the class of nonequilibrium dynamical systems which exhibits an Absorbing 
Phase Transitions(APT). This approach has better help to understand the main branches 
of the obtained fundamental diagram especially the second branch. Indeed, in this range 
of density the strong correlation between the movement of both types of vehicles depend 
on p s and ffast-, as a consequence the behaviour of the order parameter fd and the current 
J is described by the power laws(16-17-18). The critical exponents j, 5, a, v, z, rj are not 
universal because they depend on the parameters of the systems. However, the model 
shows a universel scaling law 5 = v . Furthermore, we have extended the model to the 
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case of three kinds of vehicles in which we have shown that the fundamental diagram 
exhibits three maximal currents; such result can be considered as an open question for 
experiments. 

Secondly, we recall that our aim is to propose a CA traffic flow model which we hope 
can contribute to a better understanding of the traffic flow. Hence, according to the 
previous discussions in section 3 we can say that our proposed model can give a pertinent 
indications about empirical results. Indeed, Figs.l show the existence of two branches 
with positive slops which can coexist at the same density interval as shown in Fig.8. Such 
results was obtained in empirical work [3]. moreover, in real traffic the fractions of each 
types of vehicles is unpredictable. Thus, for a fixed density C, there is an infinity of 
couples (f fast, f slow) which can give different values of the current. Thus, according to 
this observation we have studied the case of three types of vehicles where their fractions 
(f fasti, ffast2, fslow) are chosen randomly. In this situation we have combined the disorder 
in the maximal speed, the disorder in the fraction of cars and the condition under which 
we can have the metastability. Hence, the results displayed in Fig.9 look similar to the 
empirical data [7, 8]. 
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Figure captions 



Fig.l-a: The fundamental diagram (J,C) for two types of vehicles in deterministic case 
(p s = 1) with (Vmaxl = 5,ff ast = 0.3) for typel, and (Vmax2 = I, f s i ow = 0.7) for 
type2. 

Fig.l-b: The fundamental diagram (J,C) for two types of vehicles in deterministic case 
(p s = 1) with (Vmaxl = 5,ff ast = 0.6) for typel, and (Vmax2 = I, f s i ow = 0.4) for 
type2. 

Fig.l-c: The fundamental diagram (J,C) for two types of vehicles in deterministic case 
(p s = 1) with (Vmaxl = 5,ff as t = 0.75) for typel, and (Vmax2 = l,f s i ow = 0.4) for 
type2. The inset show the agreement between the analytic formula (14) (circle dots), and 
numerical similation (up triangles). 

Fig. 2: The mean speed V versus the density C for two types of vehicles, circle dots cor- 
respond to type2 with (Vmax2 = 1, f s iow = 0.25), whereas down triangles correspond 
to typel (Vmaxl = 5,ff as t = 0.75). These curves were obtained in deterministic case 

(Ps = 1). 

Fig. 3: The fundamental diagram (J,C) for two types of vehicles with (Vmaxl = 5, ff as t = 
0.75) for typel and (Vmax2 = 1, f s i ow = 0.25) for type2; circle dots correspond to de- 
terministic case p s = 1, square dots (p s = 0.7) and up triangles (p s = 0.5) correspond to 
stochastic case. 

Fig. 4: The fundamental diagram (J,C) and the mean speed V versus the density C (in- 
set) for three types of vehicles in deterministic case (p s = 1). With (Vmaxl = W,f fasti = 
0.3) for typel(down triangles), (Vmax2 = 4:,ff ast 2 = 0.4) for type2 (circle dots), and 
(Vmax3 = 1, f s i ow = 0.3) for type3 (up triangles). 

Fig. 5: The fraction of fast vehicles which succeed to overtake versus the density C. In 
the main plot, the simulation was done for two types of vehicles with (Vmaxl = 5, ffast = 
0.75) for typel, (Vmax2 = l,f s i ow = 0.25) for type2, and for three values of p s : up tri- 
angles correspond to p s = 0.9, circle dots correspond to p s = 0.8, and p s = 0.6 for square 
dots. We remark that the analytic formulas (10- 12) are verified. 

However, in order to show the effect of ffast on the behaviour of f^ (inset), the simulation 
was done in the deterministic case p s = 1, with V(maxl = 5 for typel, Vmax2 = 1 for 
type2, and for three values of ff as t'- U P triangles correspond to f fast = 0.6, circle dots 
correspond to ffast = 0.7, and ff as t = 0.75 for square dots. 

Fig.6: The scaling plot of the coexistence curve of the order parameter fa for three values 
of ff as t- square dots correspond to ff ast = 0.75, circle dots correspond to ff ast = 0.7, 
and up triangles correspond to ff as t = 0.6. The case when ff as t < 0.6 is displayed in the 
inset. We note that the scaling is valid only in the singular part of the order parameter 
fd -i.e, for C max \ < C < C m i n . The simulation was done for two types of vehicles in the 
deterministic case p s = 1 with Vmaxl = 5 for typel, and Vmax2 = 1 for type2. 
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Fig. 7: The scaling plot of the coexistence curve of the current J for C < C m i n . In the main 
plot The simulation was done for two types of vehicles in the deterministic case p s = 1 
with Vmaxl = 5 for typel, Vmax2 = 1 for type2, and for three values of ff as t : square 
dots correspond to ff as t = 0.75, circle dots correspond to ff as t = 0.7, and up triangles 
correspond to fj as t = 0.6. In this case a = 1.631 and v = 0.773. 

However, in the inset The simulation was done for two types of vehicles with a fixed 
fraction ff ast = 0.75 and for several values of p s : square dots correspond to p s = 0.3, 
circle dots correspond to p s = 0.4, up triangles correspond to p s = 0.5 and down triangles 
correspond to p s = 0.7. In this case z = 0.302 and n = 0. The plot does not contain the 
region when C > C m i n , because the system is in the absorbing phase where fd = -i.e., 
the movement of typel is dictated by the one of type2, which leads to the well known 
fundamental diagram of one species of vehicles with V max = 1. 

Fig.8: The fundamental diagram (J,C) for two different initial configurations: Up trian- 
gles correspond to homogeneous initial configuration and down triangles correspond to the 
mega jam. The simulation was done for two types of vehicles with (Vmaxl = 5, ff as t = 
0.6) for typel, (Vmax2 = 1, f s i ow = 0.4) for type2 and p s = 0.7. 

Fig.9: The fundamental diagram (J,C) and The mean speed V versus the density C (in- 
set) for three types of vehicles in deterministic case (p s = 1) with Vmaxl = 10 for typel, 
Vmax2 = 4 for type2, and Vmax?, = 1 for type3. For each value of the density C the 
fractions (f fasti, ffast2, f slow) are chosen randomly. 
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